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Abstract 



We show that perturbation theory provides two distinct mechanisms for the 
power like growth of hadronic cross sections at high energy . One, the leading 
BFKL effect is due to the growth of the parton density, and is characterized 
by the leading BFKL exponent ui. The other mechanism is due to the infrared 
diffusion, or the long range nature of the Coulomb field of perturbatively mass- 
less gluons. When perturbative saturation effects are taken into account, the 
first mechanism is rendered ineffective but the second one persists. We suggest 
that these two distinct mechanisms are responsible for the appearance of two 
pomcrons. The density growth effects are responsible for the hard pomeron 
and manifest themselves in small systems (e.g. 7* or small size fluctuations in 
the proton wave function) where saturation effects are not important. The soft 
pomeron is the manifestation of the exponential growth of the black saturated 
regions which appear in typical hadronic systems. We point out that the non- 
linear generalization of the BFKL equation which takes into account wave func- 
tion saturation effects (" pomeron loops" ) provides a well defined perturbative 
framework for the calculation of the soft pomeron intercept. The conjecture of 
a perturbative soft pomeron is consistent with picturing the proton as a loosely 
bound system of several small black regions corresponding e.g. to constituent 
quarks of size about 0.3 fm. Phcnomenological implications of this picture are 
compatible with the main qualitative features of data on p-p scattering. 
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1 Introduction 



It has been long recognized that at asymptotically high energies hadronic cross 
sections are dominated by soft nonperturbative physics. In particular the va- 
lidity of the Froissart bound for the total cross section requires a mass gap 
in the spectrum. The corresponding generation of the pion mass in QCD, or 
the generation of the glueball mass in pure gluodynamics, is a bona fide non- 
perturbative effect. However, even though the asymptotics is expected to be 
nonperturbative, perturbative dynamics may well play an important role in the 
preasymptotic regime. 

Indeed, cross sections for small objects, like highly virtual photon 7*, or heavy 
quarkonium ("onium") are perturbative up to very high energies, as long as the 
size of the system remains small. The main perturbative mechanism that drives 
the growth of cross sections with energy is the BFKL evolution It predicts 
that cross sections grow exponentially, a oc s u , where to = 41n2iV c ^ to leading 
order in a s . This behaviour is nonunitary, as it violates the Froissart bound 
a < -^t 2 . Here t is the rapidity, t = ln-%, and m denotes the mass of the 
lightest particle in the theory. The main reason for this lack of unitarity is the 
growth of the partonic density in the hadronic system. As the partonic density 
reaches a critical value of order l/a s , the BFKL approach ceases to be valid, 
and one must take into account finite density effects. 

This has been recognized and forcefully advocated in the pioneering works Q 
and later in ||. A more complete approach to QCD evolution at finite density 
has been developed more recently in Q|. These papers treat the QCD evolution 
to leading logarithmic approximation, but to all orders in the gluon density. The 
resulting system of evolution equations is a set of functional equations, and its 
study so far has not been feasible beyond the double logarithmic approximation. 

There exists however a regime in which the general nonlinear evolution simpli- 
fies. This happens when a small object (such as highly virtual photon) scatters 
on a large target (such as a large nucleus). In this case the leading nonlinear 
corrections are due to the fact that the projectile wave function at high energy 
has a large multigluon component. As soon as the number of gluons in these 
multigluon states becomes large, one has to account for the possibility that more 
than one gluon scatters, even if the gluonic density may still be small. When 
the target is large, so that the scattering probability of a gluon is parametrically 
larger than a s (for a large nucleus of atomic number A, it is 0(A 1 ^ 3 a s )), these 
corrections become important earlier than those due to high density effects. 

The system of evolution equations which takes into account these multiple scat- 



tering effects has been derived first in 
independently in the dipole picture of 



. In |3| its large N c limit was derived 
The advantage of the large N c limit is 



that the (otherwise infinite) hierarchy of evolution equations closes and becomes 
a single equation for the gluon density (or dipole scattering probability). We 
will refer to this nonlinear equation as the BK equation. The relation between 
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the equations of and || has been discussed in B, where it was shown that 
the latter is the limit of the former when the induced field density is small. This 
result was later rederived by similar methods in jj^Q 

The range of validity of the BK equation is wider than that of the BFKL evo- 
lution in the sense that it can be applied to scattering on large targets. If 
the target is large so that the scattering probability of a given probe on it is 
of order unity, the BFKL evolution violates unitarity very quickly, at rapidi- 
ties of order tsFKL oc l/a s N c . This violation of unitarity stems from the 
scattering probability rising above one locally in impact parameter space. In 
contrast to BFKL, the BK evolution ensures that the scattering probability at a 
given impact parameter is always below unity. Nevertheless, as we will discuss 
in the following, for scattering of a (small) projectile of transverse size x$ on 
a (large) target of size Rq, the BK evolution violates unitarity for rapidities 
t > t BK oc l/a s N c \n(R /x ). The reason for this violation is that the value 
of the maximal impact parameter that contributes to scattering grows expo- 
nentially with rapidity. Thus even though the scattering probabilities remain 
unitary, the total cross section grows exponentially, a oc exp{e£}, and violates 
the Froissart bound. This is the main result of ]l0| . We can give a crude es- 
timate of the exponent e, but we are not able to perform a reliable analytical 
calculation at this point. However the numerical results of give e = 0.75w, 
which is indeed compatible with our rough estimate. Thus the BK exponent is 
smaller than the BFKL one, but not very significantly. 

The reason why at high enough energy the BK evolution ceases to be valid 
is that the projectile wave function becomes dense, and further evolution is 
affected by the high density effects in the projectile wave function. The onset 
of these corrections is at t oc N In As discussed below, the gluon density 
in the center of the projectile is large at energies at which large values of the 
impact parameter become dominant. Thus, finite density corrections in the 
central region are important. Their main direct effect is to slow down the 
growth of density in the central region, but this must also have a large effect on 
the growth of the scattering probability at large impact parameters. As we will 
discuss below, this growth at large impact parameter is due to the long range 
Coulomb (or Weizsacker- Williams) fields originating from the fluctuation of the 
colour charge in the central "black" region. The finite density effects inevitably 
reduce the magnitude of these fluctuations and therefore also the Coulomb fields 
which feed the periphery. We expect therefore that these corrections are likely to 
reduce the exponential e significantly. The proper calculation of the exponential 
e must involve the analysis of the full nonlinear equation, including the wave 
function saturation effects (J]. 

1 We note that although it is claimed verbally in fl that the equations of [Ij should be 
equivalent to those of H not only for small induced fields, but also in genera" the actual 
mathematical analysis of [^j does not justify this claim. The mathematical analysis of 
(up to notational differences) is equivalent to that of [M, where the origin of the differences 
between the results of [Ml and Ml has been discussed. 
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One aim of the present paper is to furnish more details to the derivations of 
[ p"o| . This is done in sections 3 and 4, after the short recourse to the BFKL 
evolution in section 2. The other aim is to reflect on possible implications of 
these two distinct perturbative mechanisms for the power like growth of the total 
cross section. The first mechanism - the growth of partonic density - is naturally 
associated with small systems. As long as the size of the system is small and also 
the number of partons not very large, one expects that the growth of density is 
indeed the leading mechanism. This suggests that the BFKL exponent (possibly 
modified by higher order corrections) dominates the scattering cross section of 
small systems. However, if a typical hadron is a dense partonic system rather 
than a dilute one, then the density growth mechanism is rendered irrelevant 
by saturation effects, ft is only the power like growth due to the perturbative 
expansion in the transverse plane which remains effective in this case. One 
is thus lead naturally to conjecture, that it is this latter mechanism that is 
responsible for the experimentally found " soft pomeron" behaviour of hadronic 
cross sections. 

Associating the soft pomeron with a perturbative phenomenon is certainly un- 
conventional. All models of the soft pomeron that we are aware of, try to 
explain the power growth of cross sections by a nonperturbative mechanism 
[ fl2"l |f3| . On the other hand, this power growth of the cross section implies by 
the extension of Heisenberg's argument |L4| the power like distribution of " mat- 
ter" in impact parameter space. Perturbation theory naturally provides such a 
mechanism, since it operates with massless gluons and thus long range Coulomb 
fields. Thus it could well be that the role of nonperturbative physics is limited to 
taming the powerlike perturbative growth in asymptotia, rather than to provide 
an additional mechanism for a power like growth in the pre- asymptotic regime!^ 

For this scenario to be viable, the basic saturated ("black") building blocks of 
hadrons must be themselves small in size, so that the perturbative Coulomb 
like structure of the gluon fields is still relevant at these distances. It has 
been pointed out repeatedly in the literature that the distance scale associated 
with the constituent quarks is much smaller than the scale of confining physics. 



Indeed, the models |L2|, 13 deal with scales of order .3 fcrmi. The picture we 
suggest therefore, is that of a proton containing three small "black" constituent 
quarks with the size about .3 fermi in the proton rest frame. As the proton 
is boosted, the radius of constituent quarks grows with half the soft pomeron 
power, r cx s 04 . The mechanism of this growth is perturbative and the power 
itself should be calculable from the full nonlinear QCD evolution equation Q. 
Section 5 of the present paper is devoted to a more detailed discussion of this 
scenario. 



2 We note that the models of |12| |l3| . when translated into our language try to provide 
an alternative nonperturbative mechanism for the growth of the density rather than for the 
growth in impact parameter space. Since the exponent governing this growth is supposed to 
be much smaller than the perturbative one, it is hard to see how this sort of mechanism can 
survive as the leading one when juxtaposed with BFKL. 
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2 The two exponents of the BFKL equation 



To illustrate the perturbative features of high energy reactions, we choose the 
example of Deeply Inelastic Scattering. The subprocess relevant to strong in- 
teraction physics is the scattering of the photon of virtuality Q 2 on a proton. 
In the parton model, the photon counts the number of charged partons in the 
proton. The scattering cross section is then 



where the first factor is the parton level cross section while 2Vj is the number 
of partons of a given species in the proton. The number of partons depends on 
both, the resolution scale Q 2 and the energy at which the proton wave function 
is probed, x — qi+ w h ■ Here W denotes the center of mass energy of the j*p 
system. One can define the phase space density of partons, cj>(k), in terms of 
which 

N(Q 2 ,x) = / d 2 b [ d 2 kcj)(b : k,x) , (2) 



IS Jk 2 <Q 2 

where k is the intrinsic transverse momentum of the parton, and b is the impact 
parameter at which the parton is found in the transverse plane. 

Although the photon directly counts the number of charged partons, i.e. quarks, 
at low x the proton wave function is dominated by gluons. The number of quarks 
is then directly determined by the gluon content of the wave function. In the 
following we will therefore concentrate on the gluons only. 

At low x, in the leading logarithmic approximation, the gluon density <fi is 
determined by the asymptotic solution of the BFKL equation]!]]. According to 
this solution, the distribution of gluons is 

. . , In 2 b 2 kko 
<p(b, k, x, Kq) oc cxpjwt — — }. (3) 

Here t = Inl/x is the rapidity, uj = 41n2A c a s /7r and a 2 = 14£(3)iV c a: s /7r 
with £(n) being the Riemann zeta-function, £(3) = 1.202 . . .. This distribution 
depends also on fco which characterizes the initial transverse momentum of the 
gluon which gave rise to <j> through the evolution to high t. Formula (||) is valid 
for impact parameters which are not too large, namely 

< ln6 2 fcfc < a s t . (4) 



The striking feature of eq.(j3j) is that even if at low rapidity to one starts with 
a single gluon, after evolution to high enough t, the density of gluons at small 
impact parameters becomes exponentially large. Within the diffusion radius 
In 6 oc at 1 / 2 , the overall scale of density is determined by the exponential factor 
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expjwi}. Thus using eq.(Q) in the gluonic analog of eq.(|l|) gives the total cross 
section which grows exponentially with rapidity. 



The high density at central impact parameters is not the only source of growth 
in eq.(|j|). The gluon density does decrease towards the peripheral impact pa- 
rameters, but it does so rather slowly. For a given transverse momentum fc, 
the density stays finite up to impact parameters of order ln(fe 2 fcfco) oc a s t. 
One can not establish the exact proportionality coefficient from eq.(|^), since 
the validity of this equation does not extend to such large impact parameters. 
However impact parameters of order ln(6 2 fc/co) = va s t with very small v are 
still covered by eq.(|^). For such impact parameters the density is still exponen- 
tially large, although parametrically smaller than inside the "diffusion radius". 
These peripheral impact parameters do not contribute to the leading BFKL 
exponential. When the density is integrated over the impact parameter plane 
to calculate the total number of gluons, one obtains the leading BFKL result 
N(t) — f cPb <fi(b,t) — aexptot which is dominated by the impact parameters 
within the diffusion radius. However the contribution of the peripheral region 
is by itself also exponentially increasing. Were we to exclude the central im- 
pact parameters from the integration region, we would still get an exponentially 
large contribution of the form exp{ea s t}. The exponent here is smaller than the 
leading BFKL one, and thus gives a negligible contribution to the total cross 
section within the BFKL framework. However, this exponent is nevertheless 
present, and its physical origin is quite distinct from the exponential growth of 
partonic density. 

The exponential growth of the cross section violates the Froissart bound which 
requires the total cross section to grow not faster than the second power of the 
logarithm of energy 

Vtotal < m2 — ' ( 5 ) 

The leading exponential growth of the BFKL cross section is of course unphys- 
ical at high enough energies. Taken at its face value eq.(Q) in conjunction with 
eq. (ph would mean that the probability for scattering of a strongly interacting 
probe at fixed impact parameter grows without bound at large energy. This is 
inconsistent with the fact that the probability can not exceed one, and violates 
unitarity of the scattering probability at fixed impact parameter. The reason 
is that a strongly interacting probe has finite probability to interact with more 
than one gluon at fixed b when the gluon density at fixed impact parameter and 
fixed resolution Q 2 exceeds l/a s . Thus the cross section is not proportional 
to the number of gluons anymore. Multiple scattering effects must be prop- 
erly taken into account in order to relate the gluon density with the scattering 
probability. 

A nonlinear QCD evolution equation that takes into account these multiple 
scattering effects has been derived by Balitsky ||. Its large N c version was 
obtained by Kovchegov j|] using Mueller's dipole model approach. In the next 
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two sections we will discuss how far this perturbative resummation goes beyond 
the simple BFKL framework towards restoring the unitarity of the hadronic 
cross sections. 

We start our discussion with considering the evolution in the target rest frame. 



BOOST 



z 




Figure 1: Boosting the qg-dipole (x,y) generates the higher Fock component 
\q(x) q(y) g{z)) with a gluon at transverse position z. In the large N c -\mi\t, this 
corresponds to the generation of dipoles (a;, z) and (y, z). 



3 The BK evolution in the target rest frame 



The essence of the BK evolution is the following. Suppose at some initial rapidity 
to we are interested in the scattering of a probe consisting of a qq dipole on a 
large hadronic target, for example a heavy nucleus. The scattering probability 
of a dipole with legs at transverse coordinates x and y on a hadronic target is 
N(x, y). Increasing the center of mass energy (or rapidity) amounts to boosting 
the dipole to rapidity t = t + St. Under boost, the longitudinal Coulomb 
field associated with the dipole acquires a transverse part, i.e., the dipole of 
transverse size x — y generates an extra gluonic component whose density is 
given by the equivalent gluon content of the Weizsacker- Williams field of the 
dipole, see Fig. [l| 

\q(x)q(y)) - A\q(x) q(y) g(z)) (6) 



with 



A 2 



OsNc 

2ir 2 



St- 



(x - yf 



(7) 



(x- z) 2 (x-y) 2 ' 

In the large iV c limit, the gluon is equivalent to a quark- antiquark pair. This and 
the global colour conservation of the QCD evolution implies that the singlet qqg 
state is equivalent to two dipoles with coordinates x, z and y, z respectively, see 
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Fig. [l] In the leading large N c approximation, the dipoles scatter independently 
of each other. Thus the probability for the scattering of the pair of dipoles is 



N(x, z; y, z) = N(x, z) + N(y, z) - N(x, z)N(y, z) . (8) 

The last negative term, N(x, z)N(y, z), is the probability that both dipoles scat- 
ter in the same collision. Such double scattering events should be counted once 
and not twice in the total cross section, and the last term in eq. (||) corrects 
the overcounting of N(x, z) + N(y, z). 

This leads to the following nonlinear evolution equation for the dipole scattering 
probability: 

[N(x, z) + N(y, z) - N(x, z)N(z, y) - N(x, y)] . 

The first three terms in eq. (^) are just the ones discussed above, while the last 
term is the "virtual" correction which ensures that the dipole wave function 
stays normalized throughout the evolution. 

The following two properties of the BK equation eq.(^) are important for our 
discussion. First, as stressed above, it takes into account multiple scattering 
corrections. Second, within this approach the projectile wave function evolves 
according to the linear evolution of the dipole model. Interactions between the 
dipoles in the projectile wavefunction are not taken into account. This is seen 
explicitly e.g. from the original derivation of || where the density of dipoles in 
the projectile wave function n\{x,y) satisfies the linear BFKL evolution equa- 
tion, 



d a s N t 
— n\ \x,y) = 



n 1 (x,y) = ^ I d 2 z (10) 



1 1 . (x — y) 2 
^ni(x,z) + — —nxijj.z) - — —m{x,y)\ 



(y — z) 2 (x — z) 2 [x — z) 2 (y — z) 2 

The nonlinearity in the BK evolution equation comes not from the nonlinearities 
in the evolution of the projectile wave function, but from the nonlinearity in the 
relation between the dipole density and the scattering probability. This is again 
given explicitly in A single dipole (xq, yo) at initial rapidity to develops at a 
greater rapidity t into a wave function characterized by the m-dipole densities, 
n m (xo,yo,to\xi,yi; x m ,y m ,t). If the single dipole scattering probability at 
to is 7 (a;, y), the total scattering probability is given by H 

N(x ,y ,t) = 2_,\mLil{xi,yi)] n m (x ,yo,t \xi,yi; ...;x m ,y m ) . (11) 

m 

The scattering probability 7 depends on the target, but not on the rapidity. 
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The account of multiple scatterings in the BK resummation eliminates the lead- 
ing mechanism that renders the BFKL evolution nonunitary. Since the scat- 
tering probability in the BK evolution is no longer proportional to the gluon 
density, the scattering probability at any impact parameter does not exceed 
unity. This is obvious from eq.(^). When the scattering probability N(x,y) 
reaches unity at all impact parameters, the right hand side of the evolution 
equation vanishes, and the probability stops growing. 

A number of numerical fl^||M|17|| as well as analytical p^],|l9l studies of eq.(^|) 
have been performed, and they all lead to the following consistent picture: Sup- 
pose one starts the evolution from the initial condition of small target fields 
(or N(x,y) <C 1 for all x,y). Then initially the evolution follows the BFKL 
equation, since the nonlinear term in eq. ([)]) is negligible. As the scattering 
probability approaches unity, the nonlinear term kicks in and eventually the 
growth stops as the RHS of cq.(^) vanishes for N(x,y) = 1. The larger dipoles 
(large (x— y) 2 ) saturate earlier, with the smaller dipoles following at later "time" 
t. The following simple parametrization (2(J of the scattering probability gives 
an adequate description of the evolution 

N(x,y) = l-exp{-(x-y) 2 Q 2 s (t)} (12) 

with the saturation momentum Q s (t) a growing function of rapidity. Thus at 
any given value of rapidity, all pairs of size greater than Qj l {£) are saturated. 

The exact dependence of Q s on rapidity is not known, but both, the numerical 
results [|l5) , |l7j and simple theoretical estimates , are consistent with the 
exponential growth of the form 

Q a {t) = Aexp{a fl Ai} (13) 

with A of order unity. This physical picture has been anticipated several years 
ago in |2l| . 

Does saturation of the scattering probability locally in impact parameter plane 
necessarily imply that the total cross section unitarizes and satisfies the Froissart 
bound? The answer clearly is negative. The Froissart bound states that the 
inelastic cross section for the scattering of a hadron (dipole) on a hadronic target 
can not grow faster than the square of rapidity 

* < ^i 2 , (14) 

where m is the mass of the lightest hadronic excitation. To calculate the inelastic 
cross section one has to integrate the scattering probability over the impact 
parameter. Thus, in the saturation regime 

a = 7rR 2 {t) , (15) 

where R(t) is the size of the region in the transverse plane in which the scattering 
probability for hadronic size "dipoles" is unity. This radius itself depends on 
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t. To satisfy the Froissart bound the radius R(t) should grow at most linearly 
with t. The question of unitarity is therefore the question about the rate of 
growth of the " black" region, and thus is completely separate from the question 
of saturation of the scattering probability at fixed impact parameter. 

As an aside we note that for the Deeply Inelastic Scattering the unitarity bound 
is somewhat different. In this case the projectile is a virtual photon. It does not 
have a fixed hadronic size, but rather is characterized by a distribution of dipole 
sizes. The perturbative wave function of the virtual photon is well known. An 
interesting property of this wave function is that for transverse photon it has a 
long tail on the small dipole side (r 2 « Q~ 2 ) 

$ 2 (r) oca em -^. (16) 

In such a projectile not all dipoles saturate at the same energy. The scattering 
probability thus is given by the integral over the dipole sizes. At every rapidity 
the main contribution comes from the dipoles which are saturated, that is those 
with sizes above The scattering probability for a virtual photon at high 

energy (by high we mean here such that Q s » Q) 7 is given by 

W)= / d 2 r$ 2 (r)oca em lnQ s /Q. (17) 

Jr 2 <Qs 2 

With the exponential dependence of Q s on rapidity this translates into^j 

iV(7*) oc a em a s \ns/s , (18) 

and therefore 

a D is oc a em a s nR 2 (t)t . (19) 

Thus the DIS cross section has an extra power of t relative to the cross section of 
a purely hadronic process. This extra power of t is consistent with the numerical 
results of jlTj]. The basic question of unitarization however remains the same: 
what is the dependence of R on rapidity ? 

While there is no doubt that QCD is a unitary theory, and therefore indeed 
R(t) oc t, there is no guarantee that the nonlinear BK equation eqs.(|3^),(p|) 
preserves this property. In fact simple considerations indicate the opposite. 
This is especially clear from Kovchegov's derivation || (see also [^2j) where the 
density of dipoles in the projectile wavefunction is explicitly determined by the 
BFKL equation. Saturation is the result of the multiple scattering of the dense 
dipole system, rather than the slowdown in the growth of the dipole density. 
Since the transverse size of a system in BFKL evolution grows exponentially 
with rapidity, there is little doubt that the BK evolution violates unitarity of 
the total cross section. 

3 Although this is entirely academic, we note that this expression is valid only at energies 
below s c such that a em a s lns c /s(i < 1. Above this energy higher order electromagnetic 
corrections must become important so that N(^y*) saturates. 
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Target x : 

Rest 

Frame *oost 



x x x 

X K„X 

X x 



(transverse) t ^ ► t = t Q + 5t 



Target ^=(x J ,y 1 ) 

Rest ^=(x ff y ) boostN ^=(x 2 ,y 2 ) 



Frame ' Si ( x 3 >yp 

t Q — — t = t Q + 5t 

(longitudinal) 




Figure 2: The BK evolution in the target rest frame. Under boost, the initial 
dipole (xo,yo) evolves into a wavefunction containing m-dipole configurations 
with density determined by the BFKL expression (^). The spread of these 
configurations in impact parameter space leads to a finite interaction probability 
even if the initial dipole (xo,yo) was at very large impact parameter. 



We now present a simple calculation that establishes this point: 

Consider the BK evolution as the evolution of the projectile ||,||2^]. Suppose 
at the initial energy the projectile is a colour dipole of size xo. It scatters on 
a hadronic target of size Ro. As the energy is increased, the projectile wave 
function evolves according to the BFKL equation. At rapidity t the density of 
dipoles of size x at transverse distance r from the original dipole is given by the 
BFKL expression (see for example jp3|): 

32 ln*S£ 16b 2 In 2 

*■ b > *) - -2 (^5§7? ex p^ - ln ^ - (20) 

with uj — 41n2A c a s /7r and a 2 = 14C(3)A c a s /7r and ((n) being the Riemann 
zeta-function. 

Once the density of dipoles at some impact parameter b becomes larger than 
some fixed critical number, the scattering probability at this impact parameter 
saturates. The exact value of this number depends on the target, but impor- 
tantly it does not depend on rapidity. Thus the total cross section is given by 
the square of the largest impact parameter at which the dipole density in the 
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projectile wave function is of order unity. In order to estimate this directly 
from eq. (p(i|), we choose the dipole size x in ( |20| ) as x — Qj 1 ^). Recall that 
according to eq. (p^) , the dipole of this size scatters with probability one, if it 
hits inside the radius of the target Rq (in this view of the evolution only the 
projectile wave function depends on energy, while the properties of the target 
at t are the same as at to)- Thus if at some impact parameter R(t) the den- 
sity of dipoles of size Qj 1 (^o) is unity, the scattering probability at this impact 
parameter is unity as well, see Fig. ||. Requiring the exponential in eq.(|20|) to 
vanish we obtain Q 

with 



a s N c a r _ L . ui 



Numerically we find 



2tt 2 

a s N, 



< =-[-1 + ^1 + 4-]. (22) 
, e = .87w. (23) 

ZTT 

Thus, as claimed we arrive at the exponential growth of the total cross section. 




Figure 3: Under boost, the dipole density of the BK equation evolves according 
to (pp|). If the density of dipoles of critical size Qj 1 (<o) is unity at some impact 
parameter &, then the scattering probability N(b) is unity as well. 



The exact value of e given in eqs.(|22j,|23|) should not be taken too seriously. The 
point is that the explicit form of the dipole density eq.( |2p|) was derived by a 
saddle point integration, and as such is valid only for In < a s t. Since this 
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condition is not satisfied by eq.(|2l|), we can not strictly speaking use the saddle 
point expression eq. (|2^). This ambiguity however affects only the numerical 
value of e and not the parametric dependence in eq. (pl|). The reason is that 
even beyond the saddle point approximation the density has the form 

1 lniS£ 
n(x ,x,b,t) cc _ exp{a a fcF(— ^)} . (24) 
x A a s t 

The relevant condition is F = 0, and thus the solution paramctrically must be 
the same as eq. (^T|). 

It is important to realize, that although we use the BFKL dipole density of 
eq. (|20|) , our argument does not assume that the scattering probability at R(t) 
is dominated by one pomeron exchange. The only assumption is, that paramct- 
rically the total unitarized probability is the same as the one pomeron one. We 
use the criterion of dipole density only as an indicator for the magnitude of the 
total probability. The total scattering probability is given by 

oo 

P(P) = 7m(z, r)P m {x , x, b, t) , (25) 

m— 1 

where P m (xo, x, b, t) is the probability to find in the projectile wave function m 
dipoles of size x at transverse coordinate b within the area of the target radius 
Rq and 7 m is the probability of the scattering of an m-dipole state. In fact 
one should also sum over all dipole sizes smaller than R . We have neglected 
this summation in eq. ( p5| ) thus somewhat underestimating the total probability. 
Since 7 m > 7, the probability is bounded from below as 



P{b)>'Y^2Pm(.xo,x ) b,t). (26) 

m— 1 

For dipoles of size Qj 1 , the scattering probability 7 is of order unity. Thus the 
only condition that we use is 

00 

Y J Pm{x Q ,x,b,t)=0{l) (27) 

m=l 

whenever 

00 

n(x ,x,6,i) = ^ rnP m (xo,x,b,t) = 0(1) . (28) 

m=l 

The only way this condition can be violated, is if the wave function is domi- 
nated (with exponential accuracy !) by the trivial configuration with no dipoles, 
even when the average dipole number is one. Although the dipole model wave 
function is known to have relatively large fluctuations, there is nothing in its 
known properties 0, [ll] to suggest such an extreme behaviour. In fact for the 
explicit exponential model used in H, ni| our condition clearly holds. 
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To summarize, in the target rest frame, the violation of unitarity by the BK 
evolution can be understood as follows: Start with a single dipole scattering on 
the hadronic target of transverse size Rq. With increasing energy the projectile 
dipole emits additional dipoles strictly according to the BFKL evolution. The 
density as well as the transverse size of the projectile state thus grows. The 
increase in density leads to increasing importance of multiple scatterings which 
are properly accounted for in the BK derivation. This ensures that the scattering 
probability saturates locally. In the saturation regime, as long as the size of the 
projectile state R(t) is smaller than the target size i?o, the cross section grows 
essentially only due to surface effects, 

a = ttRq + 2itRoXo exp 

As long as et < In the cross section is practically geometrical. However 
once the energy is high enough, so that the projectile size is larger than that 
of the target, the total cross section is determined by the former and grows 
exponentially with rapidity according to eq. (|2l]). Thus at rapidities 

t>^rr^— , (30) 

a s N c x 

the BK evolution is nonunitary and can not be applied. 

This also illustrates that the applicability of the BK evolution crucially depends 
on the nature of the target. If the target is thick enough, so that the multiple 
scatterings become important before the growth of the projectile radius does, 
and if the target is wide enough, so that saturation occurs before the projectile 
radius swells beyond that of the target, then there is an intermediate regime 
in which the inelastic cross section remains practically constant and equal to 
7ri?g. Then BK applies in this intermediate regime. However, if the target is a 
nucleon, neither one of these conditions is satisfied. Thus the tainted infrared 
behaviour of the BFKL evolution of the projectile will show up right away and 
will invalidate the application of the BK equation. 



2ti 



-et 



(29) 



4 The target evolution picture 

The BK equation (|^) is valid only in the leading approximation in 1/N C . Beyond 
the leading order the evolution for a dipole cross section does not close, but 
rather is the first in the infinite hierarchy of equations. This hierarchy was 
derived in p|. It is useful to consider its interpretation in the frame where all 
the energy resides in the target. In this frame further increase in energy leads 
to growth of the target gluon fields. The evolution equation governs the change 
in the distribution of the gluon fields in the wave function of the target. In 
the particular gauge used in the largest component of the vector potential 
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is A + . In this gauge the S-matrix for scattering of a fast fundamental projectile 
on the target fields is given by the unitary eikonal factor U(x) 

U(x)=Pexp{i J dx-T a A+(x)}, (31) 

where T a are the generators of the SU (N) group in the fundamental represen- 
tation. The first in the hierarchy of evolution equations derived in H is 
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g a < i - gt w g w>-g/ %-')Tft; ) -,). 

x(N c Tr[U^{x)U(y)]-Tr[U\x)U(z)]Tr[U\z)U(y)]) . (32) 



The averaging in eq. (|32|) is taken over the ensemble of field strengths charac- 
terizing the target, i.e., over the target wave function. In the large N c limit the 
averages in eq. ( |32"| ) factorize and one recovers eq.@ with the identification 

N(x,y) = ^TT{l-UHx)U(y)). (33) 




BOOST 
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Figure 4: The BK evolution in the target evolution picture. As long as the dipolc 
leg yis in the white region, only the evolution of the scattering amplitude U (x) 
in (|35|) is non-trivial. 



A physically appealing reformulation of Balitsky's hierarchy of equations was 
given by Weigert |2f| in terms of a nonlinear stochastic process. The unitary 
scattering amplitude U evolves under the action of a stochastic source 



dU(x) 
dt 



{x-zf 



l-W(x)U(z) 



Ms 



1 



(x — z)'- 



■Tv[T a U\x)U(z)} 



(34) 



where U (x) and U (x) are the unitary matrices ( pl| ) in the fundamental and ad- 
joint representations, respectively. The white noise £ is characterised by Gaus- 
sian local correlations 



(#(*', *')#(*", z")) = - t"W - z") . 



(35) 
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This Langevin equation gives rise to an infinite number of equations for corre- 
lators of U which coincide with those derived in Q. 



Consider then the Langevin equation formulation, eq. (34). Assume that at 
the initial rapidity to the target is black within radius i?o- This means that 
for \z\ < i?o the matrix U(z) fluctuates very strongly so that it covers the 
whole group space. Let us concentrate on the point x which is initially outside 
of this black region. The matrix U (x) then is close to unity. Thus there is no 
correlation between U(x) and U(z), and the second and third terms on the right 
hand side of e q. (p4|) can be set to zero. This is the random phase approximation 
introduced in |25|] and used later in |H| . Note that this approximation does not 
linearize the evolution. Rather it corresponds to equating the nonlinear term in 
eq.(^) to unity for z in the black region. 

As the target field ensemble evolves in rapidity, the radius of the black region 
grows. As long as the point x stays outside the black region, we can approximate 
the Langevin equation by 



d . , a s N c f 2 (x - z 



dt V 2tH J lz]<R [x-zY 

Here we did not indicate explicitly colour indices, since they are inessential to 
the argument. We have also neglected the contribution to the derivative of U 
that comes from gluons originating from the sources outside the black region. 
Those contributions speed up the growth of U, and so by omitting them we can 
only underestimate the rate of growth of the radius of the black region. The 
formal solution of equation (|36"|) is, 



a s N c f f f 2 (x - z)i 



1-U(x,t) = d-^ / dr / d 2 z) (37) 

V 27T 2 J tQ J\ z \<R{r) {X-ZY 

Squaring it and averaging over the noise term gives 

((l-U(x,t)Y) = ^[ t drf d 2 Zj^-T 2 - (38) 

Z7r Jt J\z\<R{t) \ x ~ z ) 

As long as x is outside the black region and \x\ > R, we can approximate the 
integral on the right hand side by 



J\z\<R(t) \ x ~ z ) 



and eq.(38) becomes 



{{ l-U{x.t)f) = ^- 2 fdrR\r). (40) 

Z7r x Jt 

As the black region grows, eventually it will reach the point x. At this rapidity 
the matrix U{x) starts fluctuating with an amplitude of order one. Thus, when 
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R(t) = \x\, the left hand side of eq. ( pfl[ ) becomes a number of order one, which 
we call 1/e. We thus have an approximate equation for R(t), 

-R 2 (t) = ^ fdrR^r), (41) 
or in the differential form 

|«(«) = ^.«(«). (42) 

At large rapidities therefore the radius of the black region is exponentially large 
R(t) = R(t ) exp{^e(i - t )} . (43) 



We thus recover the result of the previous section. 



The approximations leading to eq. (42) are not strictly speaking valid when 



the point x is on the boundary of the black region. First, eq. (39) is an 
underestimate of the integral, since the inequality |x| ^> R no longer holds. 
However, this approximation can lead only to an underestimate of the rate of 
growth of R. Second, not for all points z in the black region the term U(x)U' (z) 
in eq.(|34|) can be dropped. This however is also unimportant, since when x is 
on the boundary of the black region, although the factor (1 — U(x)W(z)) is 
not striclty unity, it is still of order one for all points z. It is in fact different 
from unity only for points z in the vicinity of x. Thus, although we can not 
determine the exact numerical value of e, the functional form of the solution as 
well as its parametric dependence is given correctly by eq. (^3|). 

From the point of view of the evolution of target fields the violation of the 
unitarity can be interpreted as follows. The RHS of eq. (Q) is nothing but the 
total Coulomb (Weizsacker- Williams) field at point x due to the colour charge 
sources at points z. The stochastic noise tells us that these colour sources are 
uncorrelated both in the transverse plane and in rapidity. For such random 
sources the square of the total colour charge is proportional to the area, and 
this is precisely the factor R 2 in eq. ([io|). The incoming dipole thus scatters 
on the Coulomb field created by a large incoherent colour charge. Because the 
Coulomb field is long range, the whole bulk of the region populated by the 
sources contributes to the evolution and leads to rapid growth of R. If the 
field created by the sources was screened by some mass, the evolution would be 
perfectly unitary. To illustrate this point, let us substitute the Coulomb field 
(x — z){/(x— z) 2 in eq. (|m|) by an exponentially decaying field mexp{— m|ir— z|}. 
It is straightforward to perform now the same analysis as before. Eq. ( p9[ ) is 
now replaced by 

/d 2 zm 2 exp{— m\x — z\} = exp{— m\x — R\} . (44) 
,z\<R(r) 
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This leads to the substitution R 2 — > expjrni?} in all subsequent equations with 
the end result that 

R(t)=a a —t, (45) 
m 

which in fact saturates the Froissart bound. 

Thus the reason for the violation of unitarity is that the evolution is driven by 
the emission of the long range Coulomb field from a large number of incoherent 
colour sources in the target. 

To conclude this section we would like to discuss the relation of our results with 
numerical studies of the nonlinear QCD evolution. Studies within the framework 
of the dipole model were reported in jllj] . Ref. jll[] does not deal directly with 
the nonlinear BK equations, but rather with the onium-onium scattering in 
the framework of the dipole model. However as is clear from our discussion 
at asymptotically high energies this distinction is irrelevant. The growth of 
the transverse size of the projectile eventually determines the behaviour of the 
cross section irrespective of the nature of the target. Indeed our results are 
in agreement with those of |y]. The numerical results of [jllj clearly indicate, 
that even though the scattering probability is unitarized locally in the impact 
parameter space, the total cross section keeps on growing exponentially with 
t (Figs. 9 and 10 of Q). From Fig. 10 of Q we conclude that the power 
of the exponential is about .75w, where u> is the leading BFKL exponential. 
Interestingly, our rough estimate ( |23| ) is in reasonable qualitative agreement 
with this numerical result. 

Our results eqs. (|2l| , |42|) are in apparent contradiction with the conclusions of 
numerical work ||15| , |l6|] . The origin of this discrepancy is that these references 
solve eq. (||) within the local approximation, assuming that important con- 
tributions come only from the dipole sizes which are smaller than the impact 
parameter. Within this approximation the dependence on the impact parame- 
ter in cq. (^) becomes parametric, and the growth of the total cross section is 
determined entirely by the shape of the initial condition. The Froissart bound 
is then saturated for the exponential initial profile of N(b). The physics here 
is simple. In the local approximation of Jl5|, [ll| the gluon density of the target 
(and the scattering probability N) evolve at all impact parameters b accord- 
ing to the same exact translationally invariant equation. The density locally 
grows at all impact parameters at the same rate. This rate is the same as the 
growth of the saturation momentum and is power like with rapidity, eq. (|l3|) . 
Thus if one starts from initial configuration with the exponential density profile 
g(b) = exp{— b/ Rq}, after the evolution to rapidity t it becomes 

g(b,t) =exp{-b/Ro + a s \t} . (46) 

The scattering probability on such a system is unity at impact parameters for 
which g(b) > 1. Thus the highest impact parameter that contributes to the 
total cross section at rapidity t is b max = B ^t, and the cross section^ is a = 

4 In this discussion we neglect the dependence of N on the size of the dipole. Strictly 
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^max K t 2 ■ For initial Gaussian distribution on the other hand, the same 
argument leads to the linear growth of the cross section with t. 

This feature, namely that the asymptotic form of a is determined by the initial 
distribution is clearly an artifact of the local approximation. The reason the 
local approximation leads to this behaviour, is that it neglects the effects of far 
away black regions (where N = 1) on the scattering probability in the grey areas 
(where N < 1). As is apparent from our analysis in eq. (|3rj - ^p|) , it is precisely 
the effect of the far away black regions that drives the growth of the total cross 
section. This is due to the long range Coulomb fields originating in the central 
black region. In fact, the only contributions we kept on the RHS of eq.(|36|) are 
due to dipoles with sizes of the order of the impact parameter. In this respect 
our discussion is orthogonal to that of Jl5|, [Tg) . It is clear from the comparison 
of our results to those of jl5|, [l6) , that the effect of these long range fields on the 
total cross section is far greater than that of the local translationally invariant 
part of the evolution. Even if one starts from an exponential density profile, the 
full BK evolution generates power like and not exponential tails in the density 
at large t. These power like tails dominate the total cross section and lead to 
the universal exponential growth of a with rapidity. The local approximation is 
adequate for studying the behaviour of Q s {t) in the dense central region, as this 
is determined by local effects. It is however not a good approximation for the 
total cross section, which is dominated by the evolution of long range Coulomb 
fields. 



5 One Pomeron, two Pomeron; hard Pomeron, 
soft Pomeron 

The linear evolution of the projectile wave function inherent in the BK evolution 
is a direct consequence of the large N c limit. In this limit one can neglect the 
interactions between the dipoles in the projectile wave function. Individual 
Feynman diagrams which contain dipole-dipole interactions are suppressed by 
powers of 1/N C . Thus, the dipoles in the projectile wave function do not interact 
during the evolution and they scatter independently of each other. However, the 
number of interacting diagrams grows very fast with the number of interacting 
dipoles. As the number of dipoles in the wave function which can interact with 
each other becomes O(N), the number of the suppressed diagrams becomes 
0(N 2 ) and the suppression disappears. At high enough rapidity, where the 
dipole-dipole interactions are important, the evolution equation (^|) breaks down 
and wave function saturation effects start to play an essential role. Although 
equation (|3|) contains some 1/N C corrections, those are only "group theoretical" 
corrections reflecting the fact that at finite N c a gluon is not strictly equivalent 

speaking such dependence is of course present, and it determines the rapidity at which the 
asymptotic behaviour of a sets in. The asymptotic form of a however is independent on the 
dipole size. 
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to a qq pair. It does not contain the corrections due to the interactions of gluons 
(or dipoles) in the projectile wave function. 

The interaction probability of a dipole of size r with another dipole of similar size 
in the projectile wave function is of the order a s r 2 n(r), while the probability of 
the direct interaction with the target is 7(r). The multiple scattering corrections 
are thus more important as long as a s r 2 n(r) -C 7(1"). The projectile wave 
function corrections become important as the density grows so that 

r 2 n(r) = exp{wi} = — . (47) 
a s 

Assuming that for a large target 7 is of order one, this happens for rapidities 

toe -lnl/a, ~ 1 \nl/a s . (48) 
ijj (Jy c a s ) 

The "subleading" large N c nature of these corrections is clear if one traces 
back the explicit N c dependence of (ft8|). In fact, the coupling constant a s 
occurs in n(r) only in the combination a s N c , since the BFKL evolution neglects 
dipole-dipole interactions. Thus the a s in the numerator of eq.(fl7|) is explicitly 
subleading in N c and indeed accounts for dipole-dipole interactions. 

The above illustrates that the BK resummation improves on the BFKL evolution 
when the target is dense, that is the scattering probability 7 is larger than order 
a s . On a large nucleus of atomic number A one expects 7 cx asA 1 / 3 , and thus 
in this case BK resumms all corrections in powers of a s A^/ 3 Q. However even 
in this case the validity of the BK equation remains limited. The scattering 
probability 7 is only large at central impact parameters. At peripheral impact 
parameters the density in the target drops to zero, and so for the peripheral 
scattering events the applicability of BK is no better than that of BFKL. In 
particular, the diffusion contributions play the same role, since they are not 
associated with high partonic densities. Since at very high energies the total 
cross section is dominated by peripheral events, one does not expect the BK 
evolution to be a valid approximation in calculating the total cross section. 
Comparing eqs.(|3(]) and (fts| ) we conclude that for large targets, at rapidities at 
which the total cross section is dominated by peripheral impact parameters, the 
wave function saturation effects are important even for small projectiles. Thus 
even within perturbation theory, the wave function saturation effects are bound 
to be the ones that determine the growth of the total cross section with energy. 

Could wave function saturation effects lead to a perturbative unitarization of 
the total cross section? As seen above, the violation of the Froissart bound in 
the BK approximation is due to large and incoherent colour charge fluctuations 
in the black region. If there was a mechanism to ensure strong correlations such 
that the total colour charge in a region of fixed size L is zero, then the incoming 
dipole would feel the Coulomb field only within the fixed distance L from the 
black region. Thus the new charges produced by the evolution would only 
" split off' the edges of the black region rather than from its bulk. This scenario 



20 



is equivalent to exponential decay of the field, and would lead to a unitary 
evolution, as expected for a confining theory like QCD. However, it does not 
seem plausible that this type of correlations can be enforced on the system by 
a perturbative mechanism, nor that perturbation theory can generate any kind 
of " mass" for the gluons which would lead to exponential decay of gluon fields 
generated by sources in the black region. We thus expect that any perturbative 
corrections will generate power like tails of gluon density at large distances from 
the black region. The growth of these power tails with rapidity will inevitably 
lead to an exponential growth of the cross section. 

We expect, however, that wave function saturation effects diminish the exponent 
e significantly relative to the BK value e = .75ui. The evolution equations 
derived in Q go beyond the BK equation by including wave function saturation 
effects. They provide a well-defined perturbative framework for calculating this 
effect. 

We have seen in the previous sections that QCD perturbation theory predicts the 
existence of two distinct physical mechanisms for the power growth of hadronic 
cross sections with energy. One is due to the fast growth of the partonic den- 
sities, while the other is due to expansion of hadronic states in the transverse 
plane. The first mechanism is the leading one for systems containing a small 
number of partons. When the partonic density reaches the critical value of 
order l/a s , further growth of the density is cut off by the perturbative satura- 
tion effects. The transverse expansion mechanism is however likely to survive 
in this situation, and thus it should become the leading driving force for the 
perturbative growth of cross sections in dense systems. 

The first mechanism with the perturbative BFKL exponent of around .3 — .4 fits 
very nicely with the so called hard Pomeron utilized to fit the hadronic cross 
sections of small systems |26| . The exponent due to the transverse expansion 
has not been calculated yet, although the equations which determine it have 
been derived in ||. Although we do not know yet the value of this exponent, 
we know for certain, that it is smaller than the BFKL one. One is naturally 
lead to ask whether this perturbative growth of the transverse size of saturated 
systems underlies the experimentally observed power growth of the cross sections 
in purely hadronic processes, like pp and pp. This type of expansion, and with 
it the power like growth of the cross section should cease at asymptotically high 
energies, where the perturbation theory must become invalid. 

This picture suggests that nonperturbative effects are called for only in order 
to unitarize the cross section at asymptotically high energies, but not in order 
to furnish the mechanism for its fast growth in the pre-asymptotic regime. This 
is also natural from the following perspective. The Froissart-like behaviour is 
associated with the existence of a gap in the spectrum, as can be illustrated by 
a simple and intuitive argument due to Heisenberg jL4| . In a theory with a mass 
gap, the profile of the distribution of matter density in any target must decay 
exponentially at the periphery, p(b) cx exp{— mb}. As this target is struck 
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by a projectile, in order to produce an inelastic scattering event at least one 
particle must be produced. Assuming that the scattering is local in the impact 
parameter plane, the region of the overlap of the probe and the target must 
therefore contain energy at least equal to the mass of the lightest particle, m. 
For scattering at energy E = s/m in the frame where all the energy resides 
in the target, the target energy density is Ep(b). Thus the scattering can only 
take place for impact parameters smaller than those that satisfy E exp{— mb} = 
m. Thus b max = ^ m { s / m2 }i which is equivalent to the Froissart bound. 
Conversely, if the cross section grows as a power of energy, then the density 
distribution in the target is not exponential but power like. With p(b) oc 6~ A one 
obtains b max oc s~ . Since the power growth of hadronic cross sections persists 
in a large interval of energies, one expects that for a range of impact parameters 
the density distribution in hadronic states is power like. Perturbation theory 
provides a natural explanation of such a power like distribution. The tails of 
perturbative distribution are due to massless gluon fields emitted from the colour 
charges in the target. Even though the target is neutral, it always possesses a 
multipole moment of some order, and thus perturbatively is always accompanied 
by a long range power like tail of massless gluon field. The precondition for 




Figure 5: Schematic picture of the proton as a loosely bound system of three 
constituent quarks which provide "black" building blocks at initial rapidity to. 
Under boost, the blackness of these regions does not change since the gluon 
density is saturated, but the transverse size of the black region grows. 

applicability of this perturbative mechanism in hadronic systems is that hadrons 
themselves are built from small "black" building blocks. As explained earlier, 
the gluon fields at periphery are emitted from the bulk of the black disk, and 
not from its boundary. Thus, in order to explain the preasymptotic power like 
rise of hadronic cross sections by perturbative Coulomb like gluon fields, the 
radius of these black region must be smaller than the confining scale of QCD. 
It is in fact widely believed that QCD does naturally contain the scale of the 
right order - the scale associated with chiral symmetry breaking. In particular 
the radius of constituent quark is believed to be at most .3 fm and perhaps 
even smaller . This supports the phenomenological picture of the proton as 
a loosely bound system of three small constituent quarks. The nonperturbative 
confining force keeps these quarks confined within the spatial region of radius 
.8 - .9 fm. 
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A quantitative analysis of the non-linear evolution equations [|] which could 
substantiate our proposal of the perturbative soft pomeron, is not available 
yet. There are however some qualitative consequences of the above picture 
that can be compared to experimental data. Let us check first whether the 
value of the cross section in our model is in rough agreement with experiment. 
The total cross section of pp scattering at the lowest energy where the Reggeon 
contributions are not important, yfs oc 50 GeV is 45 mb. The total cross section 
for the scattering of two black disks of diameter d is 2nd 2 . Thus in the simplest 
model in which the proton is composed of three completely "black" constituent 
quarks of diameter d, the total pp cross section is 3 x 3 x 2nd 2 . Equating this 
to 45m6 we find d = .28 fm. This should be considered as a lower bound on 
the value of d. It is more likely that the constituent quarks are only black in 
the center and have grey peripheral regions, see Fig. ^[ For peripheral qq events 
the scattering occurs with probability / < 1. Incorporating this roughly as an 
average "greyness" factor in the formula for the cross section gives d = .28/ -1 / 2 
fm. For / = .5 we have d = .4 fm. We may thus think of the proton as a 
collection of three loosely bound constituent quarks each described by a disk of 
d = .3 — .4 fm which is essentially black in its center but grey at its boundary, 
see Fig. ||. 

This picture is quite remarkable, since the "active" area inside the proton in our 
model is much smaller than the proton radius and thus one may have worried 
that the model will underestimate the total cross section. This however does 
not happen. Further support for this picture comes from the ratio which 
is very close to 2/3. The pomeron contributions to the total hadronic cross 
sections are parametrized as a pp = 21.70s 808 and a„ p = 13.63s ' 808 @. This 
is consistent with tt having two constituent quarks. 

Another feature of the model is that although the quarks are black, the proton 
itself is not. Thus we expect the ratio of the elastic to total cross section to be 
well below the black disk value of 1/2. The experimental value of this ratio at 
t/s = 50 Gev is indeed just below 1/5 [ p8j . 

Another global characteristics of the scattering is the radius of the proton as 
measured via the shrinkage of the elastic peak, oc exp{-^-i}. For the elastic 
pp scattering this gives the value of R consistent with the proton radius R = .8 
fm | p8[ . Again at first sight this sounds like trouble, since it is much larger than 
the radius of a constituent quark. However one should realize that the radius 
of the individual quark is not relevant for this particular quantity. The process 
is elastic if the proton as a whole emerges from the interaction intact. Thus 
indeed it is the radius of the proton and not of a quark that should determine 
the i-dependence of the elastic cross section in our model. The radius of the 
quark may emerge in a similar way in the processes "elastic" with respect to 
individual constituent quark scattering. Whether one can define a subset of 
final states that correspond to such a process is an interesting question, but at 
present it is not clear to us how to do it. 
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The picture of the proton as built from three small constituent quark was en- 
voked in the nonperturbative model for the high energy scattering in The 
physics of |T3] is however quite different from that of our proposal. The quarks 
themselves in are not thought of as being black, and the growth of the 
quark-quark cross section is due to the increase of the density of the gluon cloud 
surrounding an individual quark, rather than to the increase in its transverse 
size. Ref. |D| also appeals to the scale of .3 fermi although not explicitly in 
connection with the size of constituent quarks. Again, however the mechanism 
of the growth of the cross section in Q appears to be the same as the leading 
BFKL mechanism, that is the growth of gluonic density. The role of instanton 
effects in |HJ is to limit the gluon emissions only to within the transverse sizes 
smaller than .3 fermi, and thus to cut off the expansion in the transverse plane. 
In contrast our picture assumes that nonperturbative effects (perhaps instan- 
tons) are responsible for the buildup of a black gluon cloud around each quark 
at low energy. The subsequent evolution in energy is dominated by perturbative 
swelling of these black regions. Eventually, when the size of the quarks reaches 
confining scale, other nonperturbative effects kick in and cut off further power 
like growth of the cross section. The physics of these nonperturbative effects is 
presumably the physics of confinement. 

We note also that the two distinct physical mechanisms for two pomerons have 
a consequence that they will appear with different probability in different pro- 
cesses. For example in pp scattering there are rare proton configurations which 
do not have the typical hadronic structure, but contain quarks closely bunched 
together in coordinate space. Those are the configurations responsible for the 
colour transparency effects in the hadronic scattering j2^| . These configurations 
are rare and short lived. Thus they do not have time to develop the gluon 
clouds around the quarks that could make them " black" . These configurations 
will predominantly evolve towards the increase in density, and thus will have 
energy dependence of the hard pomeron. This is consistent with the findings of 
[ ^6| that hard pomeron is present already in purely hadronic processes. 

On the other hand one does not necessarily expect the soft pomeron to appear 
in DIS even at low Q 2 . The reason is that even though at low Q 2 the DIS cross 
section has large contributions from the photon fluctuations into the states of 
hadronic size, these states do not live long enough on hadronic time scale, and 
thus do not have time to develop dense gluon clouds around the quarks. The 
energy dependence of such large but dilute states would then be of the hard 
pomeron nature. This appears to be consistent with the recent results which do 
not require the soft pomeron to fit DIS data even at very small Q 2 [ |30| . 

Perhaps the most appealing feature of this scenario is that it gives hope to 
understand the soft pomeron within the well defined, bona fide perturbative 
framework. The equations that resum the wave function saturation effects has 
been derived in [Q. Even though their numerical study is probably much more 
involved than that of cqs.(|9||32]), we think the question is interesting enough to 
motivate such an undertaking. 
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